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ABSTRACT 

We study the chemodynamical evolution of elliptical galaxies and their X-ray 
and optical properties using high-resolution cosmological simulations. Our Tree N- 
body/SPH code includes a self-consistent treatment of radiative coohng, star for- 
mation, supernovae feedback, and chemical enrichment. We present a series of ACDM 
cosmological simulations which trace the spatial and temporal evolution of abundances 
of heavy element in both the stellar and gas components of galaxies. A giant elliptical 
galaxy formed in the simulations is quantitatively compared with the observational 
data in both the X-ray and optical regime. X-ray spectra of the hot gas are con- 
structed via the use of the vmekal plasma model, and analysed using XSPEC with 
the XMM EPN response function. Optical properties are derived by the population 
synthesis for the stellar component. We find that radiative cooling is important to 
interpret the observed X-ray luminosity, temperature, and metallicity of the hot gas 
of elliptical galaxies. However, this cooled gas also leads to excessive star formation at 
low redshift, and therefore results in underlying galactic stellar populations which are 
too blue with respect to observations. Time variation and radial dependence of X-ray 
properties and abundance ratios, such as [0/Fc] and [Si/Fc], of the X-ray emitting 
hot gas are also discussed. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: 
evolution — galaxies: stellar content 



1 INTRODUCTION 

Historically, galactic astronomy was based on optical ob- 
servation. At optical wavelengths, general properties espe- 
cially for nearby elliptical galaxies are well studied. For 
example, the various scaling relations for nearby elliptical 
gala xies, such as th e Colour-Magnitude Relation ( CMR : 
e.g. 'Sandage ' 'l97^; 'Aaronson. Persson . fc Frogel I Il98lt 
[ Bower. Lucev. fc Ellis 1992a b) and fundamental p lane (e.g. 
loiorgovski fc Davis I Il987rll)ressler et al. I Il987ti . provide 
strong constraints on the formation paradigm of elliptical 
galaxies. On the other hand, progressive developments of 
telescopes and instruments have opened up other wave- 
lengths, such as X-ray and radio, and made it possible to 
observe galaxies at those wavelengths with similar or even 
higher sensitivity and resolution than those in the optical 
observations. As a result, those multi-wavelength observa- 
tions throw a new light on the study of the formation and 
evolution of galaxies. 

As for elliptical galaxies, the discovery of an X-ray emit- 
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ting hot interstellar medium (ISM) has dramatically re- 
newed the view of their physical properties and formation 
histo ry (e.g. lForman et al. ll979l : fForman. Jones, fc Tucker I 
^8^. The optical observation represents mainly the prop- 
erties of stellar component. On the other hand. X-ray ob- 
servations offer a unique tool for the study of the hot 
gas compon ent. Since there is little cold gas in elliptical 
galaxies (e.g.lBregman. Hogg, fc Roberts ll99j:lEdge booJ : 
IGeorgakakis et al. IbOOlF " those two are the major compo- 
nents of elliptical galaxies. Therefore, any successful scenario 
of elliptical galaxy formation must explain both X-ray and 
optical observation. 

The aim of this paper is to present our first attempt 
to construct a successful self-consistent cosmological simu- 
lation of elliptical galaxies whose properties are consistent 
with both X-ray and optical observations. To this end, we 
carry out a series of high-resolution cosmological simulations 
using our galactic chem odynamical evolution code (GCD+: 
iKawata fc Gibson ll20o3) . The code is a three-dimensional 
tree N-body/smoothed particle hydrodynamics (SPH) code 
which incorporates self-gravity, hydrodynamics, radiative 
cooling, star formation, supernovae (SNe) feedback, and 
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metal enrichment. GCD+ takes account of the chemical en- 
richment by both SNe II and SNe la, and mass-loss from 
intermediate mass (IM) stars, and follows the evolution of 
the abundances of several chemical elements in both gas 
and stellar components. Based on the plasma model for hot 
gas and single stellar population (SSP) for stars, this self- 
consistent chemodynamical treatment allows us to derive the 
X-ray and optical spectral energy distribution (SED) for the 
simulation end-products with minimal assumptions. Taking 
advantage of this, we compare these results for an elliptical 
galaxy in our simulation with the observational data at X- 
ray and optical wavelength directly and quantitatively, and 
discuss the physical processes of the formation and evolution 
of elliptical galaxies. 

In the X-ray observation, we focus on the well-studied 
observational values of luminosity, Lx, temperature, Tx, 
and iron abundance, [Fe/H]x. In observed clusters of galax- 
ies, there is wel l-known Lx — Tx relation, Lx oc Tx^ '- 
6 = 2.5 - 3 (e.g. lEdee fc Stewart Ill99ll: IXue fc Wu 11200(1) . 
On the other hand, the simple spherical top-hat density per- 
turbation theory predicts that the virial temperature and 
density of a cluster with mass M f ollow T oc M'^/'^il + z) and 
poc (l-f-2:)^ at a redshift of z fe.g. lPadmanabhan Ill99al . At 
typical cluster temperatures (T > 2 keV), free- free emission 
dominates the cooling process, and the X-ray luminosity is 
expected to be Lx oc MpT^^^. If all clusters followed this 
simple scaling low, we would expect Lx oc at a fixed red- 
shift, which is a significantly shallower slope than that ob- 
served. Recently, the Lx — Tx relation for groups of galax- 
ies has also been examined, and is revealed to have much 
steeper slope, h ^ 5, than that for clusters of galaxies (e.g. 
iHelsdon fc Ponman l26ool : IXue fc Wu l2000l) . Thus, in lower 
mass systems, the agreement between the simple theoretical 
prediction and the observations becomes even worse, i.e., the 
X-ray luminosity of low-temperature system is unexpectedly 
small. iPonman. Cannon, fc Navarro I 1^9^ calculated the 
specific entropy of the gas for observed groups and clusters, 
and found an "entropy fioor" of S = T/nJ'* ~ lOO/i"^/^ 
keV cm~^, where Ue is the electron immber density. This 
entropy floor reduces the X-ray luminosity by lowering the 
gas density progressively more in smaller system, and hence 
reproduces the steep Lx — Tx relations seen in groups. 

One approach to explain this entropy floor 

is "pre-he ating" of the gas befo re the system 

coUapses jEvrard fc Henrv I Il99ll: iKaiser I Il99lll . 
iLlovd-Davies. Ponmem^^Camion I i200o[l estimated 
that a pre-heating temperature of ~0.3 keV is required 
to explain the entropy floor, if the energy is injected well 
before cluster collapse. If energy injection takes place inter- 
nally, i.e. wi thin collapsed halos, more energy is required , 
^1-3 keV |Loewenstein ll2000l:IWu. Fabian. Nulsen Il2000t 
iBower et al. II2001I) . Although the source of the pre-heating 
is still unclear, this scenario is also supported by numer- 
ical simulations which take int o account the dynamical 
evolution of the system (e. g. iNavarro. Frenk. fc Whitel 
Il995l: iMuanwong et al. I l2002h . An alternative mechanism 
to ex plain the entropy floor is radiative cooling. iBrvan I 
i200Q') pointed out that smaller systems, such as groups, 
have converted more of their gas into galaxies than larger 
systems, and the gas that goes to form the galaxies was 
lower in entropy preferentially, t hus raising the mean 
entropy of the gas that remains. IVoit fc Brvan I (1200 ll) 



further argued that cooling removes low-entropy gas during 
the formation of cluster galaxies, and showed that this 
mechanism naturally explains an entropy floor of ~ 100 
keV cm^^. Hydrodynamic cosmologic al simulations also 
support these ana l ytic a r gumen t s fe.g. IPe arce et al. II2OO0I; 
Muanwong et al. I 120011 120021 : iDave. Katz. fc WeinbergI 
20021: IValdarnini II2003I) . 

The size of the potential well of elliptical galaxies 
corresponds to the low mass end of the groups observed 
in the above studies. iMatsushita. Ohashi. fc Makishima I 
I2OOQI MOMOO) showed that the X-ray luminosities of el- 
liptical galaxies are well below that expected from the ex- 
trapolation of the cluster Lx — Tx relation. Therefore, the 
entropy floor seems to be still active on the scale of ellip- 
tical galaxies. However, the above theoretical work focuses 
on larger systems, such as groups and clusters. On the el- 
liptical galaxy scale, we must take into account details of 
complex radiative cooling processes, star formation, as well 
as SNe feedback at the formati on epoch and during the 
long time evolution of the system ([Mathews fc Loewenstein I 
ll986l:lHatto'rr Habe. fc Ikeuchi Il987ll . These processes were 
simplifled in the above theoretical studies, or were difficult 
to follow due to the resolution limit of numerical simula- 
tions. Thus, to examine the Lx — Tx relation for elliptical 
galaxies is a challenging problem, and this paper will shed 
light on this issue. 

X-ray observation also provides the chemical com- 
position of the hot gas, which should contain a fos- 
sil record of past history of the elliptical galaxy evolu- 
tion, because stellar mass-loss products and SNe ejecta 
are expected to be accumulated in the hot gas. The 
standard SNe and stellar mass-loss rates predict the 
metallicity of the hot g as to be several times higher 
than the solar value (e.g. 'Loewenstein fc Mathews I Il99ll : 
[Arimoto ct al. 1997; Brighonti fc Mathews 1999). How- 
ever, previous measurements of ISM of elliptical galax- 
ies with ASCA have shown that the iron abun- 



dance is less than the solar abundance dAwaki et all 
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IMatsushita, Ohashi, fc Makishima 2000). This discrepancy 



between the theoretical predic tions and observatio ns is 
called the "iron discrepancy" dArimoto et al. I Il997t) . X- 
ray iron abundances remain a controversial is sue. Recently, 
iBuote fc Fabian I Jl998h : iBuote I (119991. 1200^ claimed that 



some X-ray luminous galaxies have roughly solar abundance, 
employing a multi-temperature plasma model. However, 
even these observed abundances are still lower than the the- 
oretical prediction. Also, high-resolution spectra taken with 
Reflecting Grating Spectrometers (RGS) on board XMM- 
Newton show less than solar abundance lIXu et al. I [2002!: 
ISakeUiou et"ani2002t iTamura et ani2003h . 

To solve this iron discrepancy, a couple of scenar- 
ios have been sugge sted. Using a chemodynamical evolu- 
tion model, iBekki I lll998l) showed that the metallicity of 
the ISM (especially in the outer regions) becomes lower 
than that of the stellar component in some elliptical galax- 
ies formed by gas-rich disk-disk galaxy mergers, because 
less metal-enriched ISM is tidally stripped away and sur- 
rounds its merger remnant. Considering evolution after 
the ga laxy formed, using a multi-phase chemica l evolution 
model. iFuiita. Fukumoto. fc Okoshil (ll^96l ■ ll997^ suggested 
that X-ray emitting hot gas is hardly enriched, because the 
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enriched gas cools more easily and drops out of the hot gas 
phase. Unfortunately, there is no study following the chemo- 
dynamical evolution of elliptical galaxies from the epoch of 
their complex hierarchical formation to the present, based on 
the self-consistent cosmological simulation. Thus, our study 
would be important to test those scenarios in the context of 
the cosmological evolution. 

From the previous studies above, radiative cooling is 
apparently considered to be a key mechanism in understand- 
ing the X-ray properties, such as Lx, Tx, and [Fe/H]x, of 
elliptical galaxies. Therefore, we first examine the effects of 
radiative cooling, comparing the results of simulations with 
and without radiative cooling. Pre-heating is also an inter- 
esting mechanism, although the origin of the pre-heating is 
unclear. Hence, we also consider the effect of SNe feedback 
which is one of candidates of the source of the pre-heating 
as well as the internal heating within the collapsed halo 
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1999: iLlovd-Davies. Ponman. fc Caniioni" I2OO0I: 
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paper, we would like to clarify the effects of radiative 
cooling and SNe feedback on the X-ray properties of 
elliptical galaxies, such as Lx, Tx, and [Fe/H]x. 

As mentioned above, we must explain not only the X- 
ray but also optical properties. Although the previous stud- 
ies pointed out the importance of radiative cooling, they 
did not answer where the cooled gas has gone. It would be 
natural to expect that such cooled gas would be the source 
of successive star formation. However, the optical observa- 
tions, such as the CMR, provide a strong restriction on such 
recent star formation. Therefore, we test the CMR of the 
stellar component of the simulation end-products as well. 
We will show that this test gives a serious constraint on the 
theoretical model. 

The outline of this paper is as follows. Details of our 
code, numerical simulation model, and analysis of X-ray 
and optical properties are described in Section |21 In Section 
13 we show the results, and compare the simulation results 
with the observational data in both the X-ray and the opti- 
cal regime. Finally, we present a discussion and conclusions 
in Section 2] Throughout this paper, we adopt for the so- 
lar ab undance the "meteoric" values in I Anders fc Grevesse I 
J1989D. 



2 METHODS 

2.1 The Galactic Chemodynamical Evolution 
Code 

Our simulations were carried out using GCD+, our original 
galactic chemody namical evolution code. Details of the code 
resented in iKawatal 1^9^ and iKawata fc Gibson I 

The code is ess entially based on TreeSPH 

Hernauist fc Katzl Il989l : iKatz. Weinberg^ fc Hernq uist 



^96 ) , which combines the tree algorithm jBarne^^^it 



19861) for t he computat i on of the gravitational forces with 
the SPH iLucv I 119771: iGingold fc Monaghan I 1197711 ap- 
proach to numerical hydrodynamics. The dynamics of the 
dark matter (DM) and stars is calculated by the N-body 
scheme, and the gas component is modeled using SPH. It is 



fully Lagrangian, three-dimensional, and highly adaptive in 
space and time owing to individual smoothing lengths and 
individual time steps. Moreover, it includes self-consistently 
almost all the important physical processes in galaxy forma- 
tion, such as self-gravity, hydrodynamics, radiative cooling, 
star formation, SNe feedback, and metal enrichment. 

Radiative cooling which depend s on the metallicity of 
the gas (derived with MAPPINGSIIL lSutherland fc Dopita I 
il993il is taken into account. The cooling rate for a gas with 
solar metallicity is larger than that for gas of primordial 
composition by more than an order of magnitude. Thus, the 
cooling by metals should not be ignored in numerical sim- 
ulations of galaxy formation dKaellander fc Hultman Ill998l : 
iKav et al. ..2000) . 

We stress that GCD+ takes account of energy feedback 
and metal enrichment by both SNe I I and SNe la, and 
metal enrichment from IM stars (see iKawata fc Gibson I 
I2OO3I for details). The code calculates the event rates of 
SNe II and SNe la, and the yields of SNe II, SNe la, and 
IM stars for each sta r part icle at every time step, con- 
sidering the ISalpeterl (^5^ initial mass function (IMF: 
mass rang e of 0.1^60 Mr:^ ) and metallicity dependent s tellar 
lifetimes llKodama lll997t iKodama fc Arimoto Ill997ll . We 
assume that each massive star (> 8 Mq) explodes as a 
Type II supernova. The SNe la rates are calculated using 
the model proposed by iKobavashi. Tsuiimoto. fc Nomotol 
^2000). The yields of SNe II, SNe la , and IM stars 
are ta ken fr omlWo osley fc Weaver") l|l995|). Ilwamo to et al.1 
(^^9^ , and Iva^^ei^Hoel^rGroenewegen' I lll997ll . respec- 
tively. The mass, energy, and heavy elements ejected from 
star particles are smoothed over the neighbouring gas parti- 
cles using the SPH smoothing algorithm. For example, when 
the i-th star particle ejects the mass of A/sN,i, the increment 
of the mass of the j-th neighbour gas particle is given by 



AAfsN.j = —MsN,iW{rij/hi 



where 



mjW{rij/hi), 



(1) 



(2) 



and W{x) is an SPH kernel llKawata fc Gibson II200I) . The 
simulation follows the evolution of the abundances of several 
chemical elements (^H, *He,^^C, ^''N, ^"^O, ^°Ne, ^"Mg, ^**Si, 
^^Fe, and Z, where Z means the total metallicity). The self- 
consistent calculation of both chemical (especially including 
the metal enrichment from SNe la and IM stars, which were 
often ignored in previous chemodynamical simulations) and 
dynamical evolution by GCD+ makes it possible for us to anal- 
yse the chemical and physical properties for both the gas and 
stellar components, which is essential for this study. 



2.2 Cosmological Simulation and the Target 
Galaxy 

Using GCD+, we have carried out a series of high-resolution 
simulations within the adopted standard A-dominated 
cold dark matter (ACDM) cosmology (f2o=0.3, Ao=0.7, 
r2b=0.019/i~^, /i=Q.7, and (78=0.9). We use the multi- 
resolution technique, to achieve high-resolution in the inter- 
esting region, including the tidal forces from the large scale 
structures. The initial conditions for the following Simula- 
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Figure 1. Dark matter density map of a portion of tlic 43 Mpc (comoving) simulation volume {upper panels), and predicted /-band 
image (physical scale) of the target galaxy (lower panels), over the redshift range 2;=3.0 to z=0.b. 
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Figure 2. R band image of the target galaxy for Models B (left) and C (right) at z=0. Overlaid contours show X-ray surface brightness. 
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tions are constructe d using the public software "GRAFIC2" 
jBertschinger 1200 j) . All our simulations use isolated bound- 
ary conditions. Initially, we perform a low-resolution N-body 
simulation of a comoving 30/i~^ Mpc diameter sphere. Mean 
separation of the particles is 30/i~^/48 Mpc. The mass of 
each particle is 2.90 x IO^^Mq, and a fixed physical soften- 
ing of 18.0 kpc is applied. 

At redshift 2: = 0, we selected an 11.6 Mpc diameter 
spherical region, which contains a few galaxy-sized DM ha- 
los. We trace the particles which fall into the selected re- 
gion back to the initial conditions at 2: = 38.6, and iden- 
tify the volume which consists of those particles. Within 
this arbitrarily shaped volume, we replace the low resolu- 
tion particles with 64 times less massive particles. The initial 
density and velocities for the less massive particles are self- 
consistently calculated by GRAFIC2, taking into account 
the density fields of a lower resolution region. Finally, we re- 
run the simulation of the whole volume (30/i~^ Mpc sphere), 
including the gas dynamics, cooling, and star formation. The 
gas component is included only within the high-resolution 
region. The surrounding low-resolution region contributes to 
the high-resolution region only through gravity. The mass 
and softening length of individual gas (DM) particles in the 
high-resolution region are 5.86 X 10^ (3.95 x 10*) Mq and 
2.27 (4.29) kpc, respectively. 

At 2 = 0, using the friends-of-friends technique, we 
found 5 stellar systems which consist of more than 2000 star 
particles. The largest stellar system has an elliptical galaxy 
like morphology, and the kinematics shows that the system 
is supported mainly by velocity dispersion rather than ro- 
tation (the ratio of the max rotational velocity to the max 
velocity dispersion, V/a ^ 0.7). We chose this system as the 
target galaxy for this study. The rest of the paper focuses 
on this target galaxy. The total virial mass of this galaxy is 
Mvir ~ 2 X 10" Mq, similar in size to that of NGC 4472, 
a bright elliptical galaxy in the Virgo Cluster. The virial 
mass is defined as the mass within the virial radius, rvir, 
which is the radius of a sphere containing a mean density of 
178f2n *^ times the critical values, pcrit ~ 3Ho/87rG, follow- 
ing Eke. Navarro, & Frcnk (T99Si)- 

As described in Section we are interested in the ef- 
fect of radiative cooling and SNe feedback. We carry out a 
series of simulations with this initial condition, using the fol- 
lowing three models. Model A is an adiabatic (i.e. no cool- 
ing) model; Model B includes cooling and weak feedback; 
Model C mimics Model B, but incorporates stronger feed- 
back. Model B assumes that each SN yields the energy of 
10^° ergs, but Model C adopts the energy of 10*^ ergs. The 
energy in Model C is unrealistically higher than the canoni- 
cal value (10^^ ergs). We use this model as an extreme com- 
parison case, to see the dependence of the strength of the 
SNe feedback more clearly. Also note that we assume that 
all kinetic energy of SNe is turned into thermal energy on a 
scale less than our resolution limit, and thus only the ther- 
mal energy of SNe is available to affect the surrounding gas 
of our resolution. 

Fig. Q shows the morphological evolution of DM in a 
central portion of the simulation volume, and the evolution 
of the stellar component in a 200 kpc region centred on the 
target galaxy for Model B. The lower panels correspond to 
the 7 band (in observed frame) image of the target galaxy, 
which is made using population synthesis techniques as ex- 
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Figure 3. Synthesised X-ray Spectrum (error-bars) and vmekal 
model fit (solid line). 

plained in the next section. We take into account the K- 
correction, but do not consider any dust absorption. The 
galaxy forms through conventional hierarchical clustering 
between redshifts 2=3 and 2=1; the morphology has not 
changed dramatically since 2=1. Fig.|21displays the final (at 
z = Q) R band image (gray scale image) and X-ray surface 
brightness profile (contours) for Models B and C. The target 
galaxy is relatively isolated, with only a few low-mass satel- 
lites remaining at 2=0. Model parameters and final basic 
properties for each model are summarised in Tabled 

2.3 Analysis 

For all the models, we examine both the resulting X-ray 
and optical properties, comparing them quantitatively with 
observation. The following sections describe how we derived 
synthesised spectra in both X-ray and optical bands from 
the simulation results. 

2.3.1 X-ray Synthesised Spectrum and Fitting 

To derive the X-ray spectrum, we need to know the 
density, temperature, and abundances of heavy elements 
for the gas component. Due to our chemodynamical 
evolution code, the gas particles in our simulations carry 
with them knowledge of those physical and chemical 
properties. Using the XSPEC ver.11.1.0 vmekal plasma 
model jMewe, Gronenschild. fc van den Oor d] ^8E ; 
Mewe^emfii^^^an den Oord 1986; Kaas tra I 1992; 
Liedahl Osterheld. fc Goldstein Ill995i) . we derive the X-ray 
spectra for each gas particle, and synthesise them within 
the assumed apertures. We next generate "fake" spectra 
with the response function of the XMM EPN detector, 
assuming an exposure time (40 ks) and target galaxy 
distance (17 Mpc). Finally, our XSPEC fitting of this fake 
spectrum provides the X-ray weighted temperatures and 
abundances of various elements. In this section, we describe 
the process of this analysis in detail. 

First, we generate the X-ray spectrum for each gas 
particle. We use the subroutine of the XSPEC vmekal 
model. Input parameters of the vmekal model are the hy- 
drogen number density, temperature, and abundances of 
He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe, and 
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Table 1. Model Parameters. 





energy per SN 


'^vir 


Mass 


within rvir 


(M0) 


Number of Particles 


Name 


(ergs) 


(kpc) 


Gas 


DM 


Star 


Gas 


DM Star 


pa 




821 


2.14 X 10^2 


1.60 X IQl^ 





36457 


40492 


B 


1050 


821 


1.45 X 10^2 


1.60 X IQi^ 


1.08 X 10^2 


21136 


40454 18966 


C 


10^2 


825 


1.83 X 10^2 


1.62 X IQi^ 


6.02 X 10" 


30765 


41036 10782 



no cooling or star formation is involved. 



Ni, with respect to the hydrogen abundance. The temper- 
ature for each particle are given directly from the sim- 
ulation output. The hydrogen number density is derived 
by the gas density times the hydrogen abundance divided 
by the mass of hydrogen atom, i.e. the proton mass. Due 
to the limitation of the size of memory, our numerical 
simulation follows the abundance evolution of only '^H, 



^He,i2c, "N, i«0, 2°Ne, ^^Mg, ^Ssi, ="^Fe. Simply, we as- 
sume the abundance of these elements with respect to 
"^H is the same as the total abundance of all isotopes 
with respect to the total (^H and '^H) hydrogen abun- 
dance, e.g. i8o/iH=0(=i'^ O -h'^ O -H^^ 0)/H(=i H +2 H). 
For the elements whose abundances are not calculated 
in our simulation, we assume [Na/H] = [A1/H] = [0/H], 
[S/H] = [Ar/H] = [Ca/H] = [Si/H], and [Ni/H] = [Fe/H], here 
[X/H] = logio(X/H) - logio(X0/H0). This assumption 
should not be an issue, because we use the same assumption 
when we fit the synthesised spectrum, as described later. Fi- 
nally, we generate the spectrum for each gas particle within 
the energy range between 0.1 and 20 keV with 1800 bins, 
using the vmekal model. 

To compare our results with the X-ray properties ob- 
served within a certain radius, we assume an aperture size. 
The X-ray fluxes for all the gas particles within the assumed 
apetrure radius are added up in each bin, and the synthe- 
sised spectrum is stored as a FITS file to enable XSPEC 
to read it as a model using atable command. Next, we 
make "fake" spectrum using f akeit command in XSPEC. 
In this paper, we suppose the observation by EPN detector 
on XMM-Newton, and the distance of the target galaxy is 
17 Mpc. We use epnJf20_sY9_thin.rmf for the EPN response 
file, and adopt 40 ks exposure time. Any absorption compo- 
nent or background are not taken into account. We binned 
the spectra to achieve at least 25 counts per bin, using "grp- 
pha" tool of HEASOFT ver.5.1. Fig.glshows an example of 
fake spectrum with i? = 35 kpc aperture for Model B. 

Finally, using XSPEC, the generated fake spectrum 
is fitted by vmekal model within the energy range be- 
tween 0.5 and 4 keV. As mentioned above, here we assume 
[Na/H]=[A1/H]=[0/H], [S/H] = [Ar/H]=[Ca/H]=[Si/H], and 
[Ni/H] = [Fe/H]. We fix the values of [He/H], [C/H], and 
[N/H] to be solar, i.e. zero, because those abundances do 
not affect the X-ray spectrum within the energy range which 
we use. Consequently, the free parameters for the fitting are 
the temperature and abundances of [O/H], [Ne/H], [Mg/H], 
[Si/H], and [Fe/H]. As shown later, except for Model A, two 
temperature vmekal models provide a better fit than one 
temperature models, irrespective of the aperture. In a two 
temperature model, we assume that the hot and cold com- 
ponents have the same abundances, to reduce the number 
of fitting parameters. The solid line in Fig.j^shows the best 



56t 



fit two temperature vmekal model for the fake spectrum of 
Model B with an aperture i? = 35 kpc. For Model A, we set 
the metal abundance to be zero, i.e. [X/H]= —00, because, 
as there are no stars, no metals are produced, and find that a 
one temperature model can fit the spectrum very well. Such 
spectrum fitting gives us the X-ray weighted temperature, 
and the abundances of various elements. 

Due to this analysis, we take into account the line emis- 
sions, and the energy range in fitting. This is more so- 
phisticated way than a simple Bremsstrahlung-based model, 
where X-ray luminosity is assumed to be proportional to 
^2yi/2^ which is conventionally used in numerical simu- 
latio n studies of X-ray properties o f clusters of galaxies 
(e.g. iNavarro. Frenk. fc White^ll995^ . Some previous stud- 
ies us e d mekal mo d el of X SPE C (e.g. [Mathiesen fc Eyrard I 
I2OOII: iToft et al. I |2002^ or iRavmond fc Smith I l| l9r 



model (e.g. iKatz fc White I Il995t iMuanwong et al. 1 1200: 

to derive the X-ray properties. Unfortunately, they fixed 
the metallicity to be the assumed value, because of 
the lack of information about chemical evolution, and 
they did not do any fitting of the gen erated spectrum. 
iFuiita. Fukumoto. fc Okoshi l (119961 Il997^ used meka model 
of XSPEC to derive the X-ray spectrum from the results of 
their one-dimensional chemical evolution model, and get the 
temperature and abundances by fitting the spectrum using 
XSPEC. However, they assumed the solar abundance pat- 
tern, i.e. [Fe/H]=log(Z/ZQ), [X/Fe]=0, Z is the total metal 
(elements heavier than He) abundance, and Zq is the solar 
abundance. Thus, the combination of our chemodynamical 
simulation and this analysis offers a higher level compari- 
son between theoretical models and observational data than 
that in the previous studies. The dependence of the results 
on the analysis method is briefiy discussed in Section [3.31 

In the following sections, we compare the simulation re- 
suhs with the X-ray observational data in MOMOO. MOMOO 
measured the various X-ray properties, such as luminosity, 
temperature, metallicity, within an aperture radius of four 
times the B-band half-light radius. Unfortunately neither 
the half-light radius nor the optical luminosity of our models 
agree with the observed ones, because of a large population 
of young stars in the central part of the simulated galaxy (see 
Section r3.2l for details). Since we cannot fix the aperture us- 
ing the same way as the observation due to this problem, we 
show the results with different apertures, and discuss the ef- 
fects of aperture. For this purpose, 20 kpc and 35 kpc radius 
apertures are used. The former one roughly corresponds to 
four times the B band half-light radius for the simulation 



end-products, when the light from young (ta 



Gyr) 



stars are ignored. The latter one is similar to four times the 
B band half-light radius of NGC 4472 {n, ~ 8.5 kpc). These 
apertures are applied to the projected distribution of the 
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gas particles within the virial radius of the target galaxy at 
z = Q. We use an arbitrary direction for projection, and the 
same projection is used when we derive the optical proper- 
ties. We have checked that the results we show in this paper 
do not depend on the direction of the projection. 

In this paper, we show the X-ray luminosity in the 0.5- 
10 keV pass band, because MOMOO showed the X-ray flux 
in this pass band. We focus on the X-ray emission from 
the hot ISM. The X-ray spectrum derived with the above 
analysis is an idealised spectrum from the hot ISM, be- 
cause the observed X-ray flux comes from not only the hot 
ISM but also the background and the point sources, such as 
X-ray binaries. MOMOO carefully reduced the background, 
and also estimated the flux of the hard component which 
is considered to come from low-mass X-ray binaries. We 
assume that the flux of the soft component measured in 
MOMOO would be dominated by the contribution from the 
hot ISM, and use them in the following section. Nevertheless, 
it is important for this type of study to identify the X-ray 
sources, and evaluate the X-ray flux from the hot ISM cor- 
rectly using high spatial and spectral resolution observation 
iSarazin. Irwin, fc Bregman ..20Qlll . 



2.3.2 Optical Population Synthesis 

The optical photometric properties which are examined in 
this paper are derive d using the same p rocedure as that 
used and described in iKawata ] ll2001allbl) . Here, we briefly 
describe this procedure. In our simulations, the star parti- 
cles each carry their own age and metallicity "tag", due to 
the self-consistent calculation of the chemodynamical evolu- 
tion. This enables us to generate an optical-to-near infrared 
SED for the target galaxy, when combined with our popu- 
lation synthesis code adopting the SSP of the public spec- 
trum and chemical evolut ion code, KA97 tKodama ...lOQ^ 
iKodama fc Arimoto Il99^ . It is worth noting that our sim- 
ulation assumes the same shape and mass range of IMF 
and the same lifetime of stars as those assumed in the SSPs 
provided by KA97. Here, the SED of each star particle is 
assumed to be that of an SSP, i.e. a coeval and chemically 
homogeneous assembly of stars. Since the observational data 
which our results should be compared with provide the lu- 
minosity distribution projected onto a plane, we have to 
derive the projected distribution of the SED from the three- 
dimensional distribution of star particles. Finally, we obtain 
the images projected in the direction used in the X-ray anal- 
ysis. For example. Fig. |5| shows the i?-band images (gray 
scale image) for Models B and C, where a 1000 x 1000 pixel 
mesh is chosen to span a square region of side 320 kpc, and 
the flux of each star particle is smoothed using a gaussian fll- 
ter with the filter scale of 1/4 of the softening length of the 
star particle. These images provide similar information to 
the imaging data obtained in actual observations. Thus, we 
can obtain various photometric properties from these images 
in the same way as in the analysis of observational imaging 
data. In the following analysis, we use images similar to the 
ones displayed in Fig. |21 but employing a 750 x 750 pixel 
mesh to span a squared region of side 100 kpc. The global 
photometric properties, such as the total luminosity and the 
colours, are obtained from the projected image data. In this 
paper, we focus on the colour and total magnitude of the 
target galaxy. For simplicity, we assume that the luminosity 
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Figure 4. Comparison of the simulated and observed (crosses 
with error-bars) Lx — Tx relations. The observational data are 
taken from MOMOO. The triangle/circle/square indicates the pre- 
dictions of Model A/B/C. The open (filled) symbols denote the 
values evaluated within the radius of 35 (20) kpc. The dashed 
line re presents the rel ation obtained from the adiabatic simula- 
tion of lMuanwong e~a \. |200l]). Open diamo nd shows the value 
for Model A evaluated by the same method as lMuanwong et al. I 
42001,). 



within an aperture rains of 50 kpc is the total luminosity. 
The colours are measured using the same aperture size as 
the observation which we compare with. These optical prop- 
erties are discussed in Section [3.21 



3 RESULTS 

First, we discuss the global properties of the target 
galaxy in different models. Section 13.11 examines the X- 
ray properties, and compares them with the observed X- 
ray luminosity— temperature (Lx — Tx) relation and X-ray 
luminosity— iron abundance {Lx — [Fe/H]^) relation of el- 
liptical galaxies in MOMOO. Section [3.21 discusses the optical 
properties of the target galaxy. Here, we focus on the CMR, 
which is well-established by the previous observation, and 
gives strong constraints on the stellar population of the el- 
liptical galaxy. In Section 13.31 we mention the X-ray prop- 
erties at different radii, i.e. the radial gradients. Then, the 
dependence of the derived X-ray properties on the analysis 
method is briefly argued. We also show time variation of 
X-ray properties in Section Finally, abundance ratios of 
the X-ray emitting hot gas are discussed in Section [3.51 



3.1 X-ray Properties at the Central Region 

3.1.1 Lx ^ Tx Relation 

Fig. 2] shows the predicted Lx — Tx relation for the three 
models at z =0; crosses with error-bars represent the obser- 
vational data from MOMOO. As shown in Table Inland Sec- 
tion l2.3.Tl except for Model A, the two temperature vmekal 
model is applied to fit the synthesised X-ray spectrum. Thus, 
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Table 2. Model Results of X-ray Spectrum Fitting. 

Aperture i? < 20 kpc 

Lx" Fe/H fcTi kT2 EMi/EMa^ xVm'' Ng(i?)'^ Ng,hot(/?)^ 
Name (10"'i erg s'l) (solar) (keV) (keV) 

a/ 30.0 0.0 0.45j;o;oo - - 531.9/511 1607 1606 

B 0.333 0.17lg;g3 l-52lg;J° 0.73^0.0? 0-323to.J?l 260.5/273 245 146 

C 2.87 0.2ll°;°j 1-221°;°^ 0.47^0^3 0.169to.olo 479.1/486 438 372 

Aperture i? < 35 kpc 

A-'' 47.1 0.0 0.45j;o;oo - - 478.9/543 3657 3654 

B 0.710 0.12+J5;°2 l-39+o;°7 0.67+°;q4 0.451+°,!,g3 311.7/347 542 411 

C 5.48 0.12j;g;°J l-lStoloi 0-54t°;o2 0-229to:on 549.7/572 1166 1058 



" X-ray flux (0.5-10.0 keV). 

Ratio of emission measure of the cooler and hotter component, 
degrees of freedom. 

Number of gas particles within the aperture. 
^ Number of hot (T > 10® K) gas particles within the aperture. 
/ IT model fit is applied, and abundance is fixed to 0. 



Fig.^Jplots the mean temperatures weighted by the emission 
measure of the two components for Models B and C. 

The adiabatic model (Model A) appears incompati- 
ble with the observational data due to its excessive lu- 
minosity and low temperature. The X-ray luminosity is 
too high, and/or the X-ray weighted temperature is too 
low. The dashed line of Fig. 2] shows the Lx — Tx rela- 
tion obtained from the a diabatic cosmological simulation of 
iMuanwong et al. I J200ll) . The result of Model A is rather 
close to this line than the observational data. Note that 
iMuanwong et all i200ll) derive the luminosity with a simple 
analytic formula, and do not apply any aperture or X-ray 
pass band, i.e. the summation of the cooling rates for all 
the gas particles with temperature larger than 10^ K. Thus, 
their luminosity means the bolometric luminosity, which is 
the reason why the luminosity for Model A is smaller than 
the dashed li ne. The way to der i ve tem perature is also differ- 
ent between IMuanwong et ai~l (|2Q^) and us. The temper- 
ature derived by our analysis should be higher than those 
by their analysis, because of the aperture and X-ray pass 
band which we apply. In fact, the open diamond in Fig. 0] 
shows the result for Mo del A derived by the same way as 
IMuanwong et all (|200lf) . and are in good agreement with 
the dashed line. 

On the other hand, the inclusion of radiative cooling 
leads to lower luminosities and higher temperatures. As a 
result, models with cooling (Models B and C) are in bet- 
ter agreement with the Lx — Tx relation of the observed 
elliptical galaxies. These results are q u alitatively consis- 
tent w ith t hose of iPearce et al. I ||2000|) , I Muanwong et al. I 
J2002D. andlDave. Katz. fc Weinberg I ll200dh (bu t see also 
lYoshikawa. Jing. fc Suto ll2000HLewis et al. IbOOoll . 

Fig. |S] shows the effect of cooling more clearly. The left 
panels show the density profiles of gas and DM. The DM 
profiles are the spherically averaged density profiles. We con- 
firm that the DM density profile for Model A is similar to 
the un iversal profile suggested by Navarro. Frenk. fc White 
il997l NFW profile) for this virialised mass halo and the cos- 
mology which we adopted. The concentration parameter, c, 
of NFW profile is a function of the mass of the system, 
which is represented with M200, the mass within the radius 



of a sphere containing a mean density of 200 pcrit- Using 
NFW routines provided by Julio Navarro^, we obtain the 
concentration parameters of c = 6.72,6.71,6.70 for Models 
A, B, and C, whose masses are M200 ~ 1.62, 1.66, 1.67 x 10^^ 
Mq, respectively. We also ran an N-body only (Sib = 0) 
simulation using the same initial condition, and confirmed 
that the DM density profile is identical to that in Model 
A. Thus, the inclusion of the adiabatic gas dynamics does 
not change the dynamics of the DM. On the other hand. 
Models B and C have a steeper DM density profile in the 
central region (r < 20 kpc). This is due to the gravitational 
drag by the central baryon concentration which is induced 
by radiative cooling. Thes e conclusions are consistent with 
the previous studies (e.g. iLewis et al" 1 l200d : I Pearce et al. I 
1200(1: lYoshikawa. Jing. fc Sutoll20n^^ ■ ahhough those pre- 
vious studies focus on larger systems, such as clusters of 
galaxies whose mass is Mvir > 10"Mq. 

The gas density profiles are shown with dots for all the 
gas particles within the virial radius. The gas density pro- 
files are smooth, except for some particles with high density 
at large radii which are due to gas within the satellite galax- 
ies seen in Fig. (21 In Model A, the gas density follows the 
DM density in the outer region. However, in the inner re- 
gion, the gas density profile becomes flatter than that of 
the DM, and the gas densities are smaller than the DM 
density. The radius where flattening happens (~ 20 kpc) 
is much larger than the softening length. Therefore, this is 
considered to be caused by the difference in dynamics be- 
tween DM and gas whose t hermal pressure keeps the equilib- 
rium density lower (sec also'Navarro^^^nk. ^^V hite Il995l : 
[Eke. Navarro, fc Frenk 1998; Pea rce et al. II2OOOI) . 

On the other hand, the inclusion of radiative cooling in 
Model B changes the gas density profile dramatically. The 
gas density profile no longer follows the DM density profile 
within a radius of 200 kpc. Between this radius and 10 kpc 
the gas density is much less than the expected density from 
the DM density [puM^b / (S^o — ^b)), and only recovers that 
density within 10 kpc. Consequently the central gas density 
for Model B becomes higher than that for Model A. This is 

^ http:/ /pinot.phys.uvic.ca/~jfn/charden/ 
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Figure 5. Density vs radius (left) and density vs temperature 
(riglit) distributions of gas particles for Model A (upper), Model B 
(middle), and Model C (bottom). In the left panels, the gray thick 
solid lines show the DM density profiles, and the gray dotted lines 
show the NF W profiles (see text for details) . The DM densities are 
scaled by f2i,/(r2o — f!;,) to compare with the gas density profiles. 
The solid curves in the right panels separate the region where 
the cooling time (t(.Qoi) is shorter (upper region) and longer than 
the Hubble time (tn). In the right panels, the black (gray) dots 
represent the particle within (outside of) the projected radius of 
35 kpc. The vertical gray solid lines show the mean temperatures 
weighted by emission measure of the two temperature components 
(described by dotted lines) obtained by the spectrum fitting. Since 
the spectrum for Model A is fitted by one component, only the 
solid line is appeared. 



the effect of radiative cooling, and the right panels of Fig. |^ 
show this effect clearly. The right panels of Fig. |5] show the 
density and temperature of the gas particles. The lines in 
the right panels show the density and temperature where the 
cooling time (tcooi) equals the Hubble time (tn) at z; = 0, 
i.e. the gas above this line can cool within Hubble time. In 
the adiabatic case of Model A, the hydrostatic equilibrium 
condition of the gas is determined independently of this line. 
Comparing Model B with Model A, the inclusion of radia- 
tive cooling ensures that the temperature of the gas in the 
upper region of the line decreases to T < 10* K where the 
cooling becomes inefficient. In other words, cooling changes 
hot gas to cold gas. The cold gas loses thermal energy and 
falls further into the potential well. This cold dense gas cor- 
responds to the high density gas within 10 kpc in the left 
panels for Model B. This conversion of the hot gas to the 
cold gas decreases the density of the hot gas in the outer re- 
gion, and makes less dense hot gas in the radius between 10 



and 200 kpc in Model B. Because the cold gas (T < 10"* K) 
does not produce any emission on the observed X-ray band, 
this decrease in the hot gas density is the reason why radia- 
tive cooling drives down the X-ray luminosity. Moreover, to 
keep the cooling time long, the denser gas has to be hotter 
in Model B. This is the reason why radiative cooling drives 
up the X-ray weighted temperature. 

On the other hand. Model C has similar Lx and Tx to 
Model B. These are caused by the effect of radiative cooling 
described above. However, Model C has slightly lower Lx 
and higher Tx than Model B. The reason for this is also seen 
in Fig.|S] Comparing Model C with Model B, some fraction 
of the hot gas in Model C stays above the line of tcooi = tn. 
This is due to the stronger heating by SNe which is balanced 
with radiative cooling, and makes the net cooling time for 
this gas longer. Therefore, the hot gas in Model C is allowed 
to be slightly denser and cooler than those in Model B. As a 
result, Lx (Tx) for Model C becomes slightly higher (lower) 
than those for Model B. 

Fig. |1] shows the X-ray luminosities and temperatures 
measured in the different apertures. In Model A, tempera- 
ture is independent of the aperture, which means that the 
X-ray emitting hot gas has little radial temperature gradi- 
ent within 35 kpc, as seen in Fig. |S| On the other hand, in 
Models B and C, the smaller aperture results in the higher 
temperature. Fig.j^also shows that the hot gas (T > 10^ K) 
with higher density has higher temperature. In other words, 
there is a negative temperature radial gradient in the models 
including radiative cooling. We will come back to this issue 
in Section 

Finally, we comment on the relationship between the ac- 
tual temperature distribution of the hot gas particles and the 
temperatures obtained by the spectrum fitting. Right panels 
of Fig.|S]show the two temperatures {kTi for higher temper- 
ature; kT2 for higher temperatuer) as well as the mean value 
{kTm) obtained by the spectrum fitting within the projected 
radius of i? = 35 kpc. These panels show that kTm corre- 
sponds to the temperature of the hot gas with high density, 
and fcTi and kT2 represent the range of the temperature 
of the hot high-dense gas. Although there are a significant 
amount of the gas particles with the temperature less than 
kT2, their density is too low to contribute to the X-ray spec- 
trum. Fig.|S]also shows that the two components of kT\ and 
kT2 are induced by the projection effect. In fact, we con- 
firmed that a single temperature model is good enough to fit 
the X-ray spectrum from the gas particles within the three- 
dimintional radius of r = 35 kpc ^ . This is consistent with 
the deprojected a nalysis of the recent high spatial resolu- 
tion obser vations (' Xu et al. iboOSi: iMatsushita et al. II2002I : 
iBuote et al. .■200a') . 



3.1.2 Lx - [Fe/H]^ Relation 

Fig. |S| compares the X-ray weighted iron abundance of our 
simulations with the observational data of MOMOO. As the 
adiabatic model of Model A does not form any stars (not 
having any cooling), we show the results only for Models B 



^ Throughout this paper, we use lower-case r for three- 
dimensional radius, and upper-case R for projected two- 
dimensional radius. 
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Table 3. Mass fraction of the ejected iron captured by each component, and mass fraction of each component 
within r = 35 kpc. 

Name star cold gas hot gas escape star cold gas hot gas 
B 0.749 0.015 0.0008 0.235 0.986 0.009 0.005 

C 0.694 0.023 0.002 0.281 0.962 0.011 0.027 



Iron mass of each component within r = 35 kpc. 

Iron mass ejected from stars within r = 35 kpc until 2 = 0. 
Mass of each component within r = 35 kpc. 
Total mass of baryon within r = 35 kpc. 
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Figure 6. Comparison of the simulated and observed (crosses 
with error-bars) Lx — [Fe/H]^ relations. The circle (square) 
shows the result of Model B (C); stars show the mean metal- 
licity of the stellar component for Models B and C. The open 
(filled) symbols denote the values evaluated within the radius of 
35 (20) kpc. 



and C. We measure [Fe/H]x = log (Fe/H)^ - log(Fe/H)Q 
using apertures with radii of 20 and 35 kpc. As explained in 
Section [2. 3. II this [Fe/H]x is derived by X-ray spectrum fit- 
ting, and means the X-ray emission weighted iron abundance 
of the hot gas. The star symbols plotted at the same X-ray 
luminosities as the other symbols show the mean (weighted 
by mass at z=0) stellar metallicity within each aperture. 
Note th at we use for the sol ar abundance the "meteoric" 
values in I Anders fc Grevesse 1 1^89). 

As described in Section^ the iron abundances observed 
in the X-ray satellite are significantly lower than the mean 
iron abundances of the stellar component. Since the iron 
abundances of the hot gas arc expected to be higher than 
those of stars (e.g. Arimoto ct al. 199[J, this is called "iron 
discrepancy", and one of the biggest mystery of the X-ray 
observation of the hot gas of elliptical galaxies. Our model 
results show lower [Fe/H]x of the hot gas, compared to their 
stellar [Fe/H]. Although our model results show even lower 
[Fe/H]x than the observed ones, these are consistent with 
the observed iron discrepancy. Therefore, it is interesting to 
see what has caused such a low [Fe/H]x of the hot gas, in 
order to interpret the observed iron discrepancy. 

To this end, we examine how much iron is held by 



each component with respect to the amount of iron ejected 
from the stars. Table |3 shows those values for the stars, 
cold gas, and hot gas, respectively. We measure those values 
within the three-dimensional radius of r = 35 kpc, which 
encompasses a significant fraction of the stellar compo- 
nent. We found that the deprojected [Fe/H]x within three- 
dimensional radius of r = 35 kpc is systematically higher 
than the projected one. However, because the difference is 
small (A[Fe/H[]j,f < 0.1), we conclude that the projected iron 
abundance within i? = 35 kpc is dominated by the contri- 
bution from the central region in three-dimensional space. 
Here, the "hot gas" is defined as gas with T > 10^ K, and 
"escape" means the remainder after subtracting the total 
amount of iron within r = 35 kpc from the total ejected 
iron from stars within r = 35 kpc. Although it is not true 
that all iron within this radius comes from the stars within 
this radius at a = 0, it gives us a rough idea where the 
ejected iron has gone. Table|21shows that the stellar compo- 
nent holds most of the iron. On the other hand, comparing 
the cold gas with the hot gas, the cold gas holds more iron 
than the hot gas, although their masses are similar. 

This is also demonstrated in Fig. |7| which compares 
[Fe/H]x obtained by the spectrum fitting with the distri- 
bution of [Fe/H], density, and temperatures for the gas par- 
ticles within the projected radius of i? = 35 kpc. The value 
obtained by the fitting, in fact, corresponds to the [Fe/H] of 
the hot high-dense gas. The cold gas, which does not con- 
tribute to the X-ray spectrum, has higher [Fe/H]. Since the 
stars formed from the cold gas, the mean [Fe/H] of stellar 
component is similar to that for the cold gas. 

This metal poor hot gaseous halo is explained by the 
following mechanism. Mass-loss from stars and SNe pref- 
erentially enrich the gas in the central region, where the 
cold gas is dominant, because radiative cooling is efficient 
(Fig.|^. Also, enriched hot gas can cool more easily than un- 
enriched gas, because th e cooling is more efficient in the gas 
with higher metallicity JSutherland fc Dopita lITggl) . Con- 
sequently, the cold gas holds more iron than the hot gas. 
Moreover, the high density cold gas is incorporated into fu- 
ture generations of stars, and thus a large fraction of the iron 
ejected from stars is locked into future generation of stars. 
As a result, the hot gaseous halo which emits X-ray is not 
enriched efficiently, leading to a lower [Fe/H]x for Models 
B and C. Hence, radiative cooling ensures that the hot gas 
has a low metallicity. 

We have analysed the mass fraction of the gas of a pri- 
mordial origin, i.e. the gas which is not ejected from stars, 
in the hot (T > 10^ K) and cold (T < 10® K) gas within 
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-R = 35 kpc. Here, we consider that the mass of the gas of 
a primoridal origin is Mg,p = ng x mg.ini, where ng is the 
number of gas particles within R — 35 kpc, and mg,ini is the 
initial mass of the gas particles, because the mass of each 
gas particle increases only by the mass deposition from stars 
in our modeling of the star formation Jl^^ligii). As a 
result, we found that 98 (98) % of the hot gas is a primordial 
origin, while 67 (63) % of the cold gas is a primordial origin 
in Model B (C). This also means that the gas ejected from 
stars is locked into the cold gas, and hardly stays in the hot 
ISM. In addition, because the hot gas is a primordial origin, 
i.e., the hot gas consists of the gas infalling from the outer 
region, the metallicity of the hot gas stays low. 

The aperture effect is also important in [Fe/H]x. Fig- HI 
shows [Fe/H]x within the different apertures, such as 20 and 
35 kpc. In both Models B and C, smaller aperture leads to 
higher [Fe/H]x, which means that the hot gas has a negative 
radial gradient of metallicity within 35 kpc. We will discuss 
the radial gradients of X-ray properties in Section 13.31 in 
more detail. 

Finally, we comment on the limitation of our analysis 
of the chemical composition for different gas components. 
We are not taking into account any thermal conduction or 
evaporation process, which is able to change the cold gas to 
the hot gas. Since the cold gas is metal rich, if a significant 
amount of the cold gas are transfered to the hot gas, the 
metallicity of the hot gas would be higher. Unfortunately, 
there is no good understanding of how such physical pro- 
cesses are implemented in the numerical simulation. How- 
ever, we consider that this process affects in a similar way 
to SNe feedback, which is also able to heat the cold gas. 
Therefore, our studies comparing different strength of the 
SNe feedback should provide us a hint of the total effect 
of the various heating processes. Also we ignore the mixing 
of heavy elements among the gas particles. If the mixing of 
heavy elements between the cold and hot gas are efficient, 
the metallicity of the hot gas would be higher. 

3.2 Optical Properties 

In the last section, we showed that radiative cooling helps 
to explain the observed X-ray properties, such as Lx, Tx, 
and [Fe/H]x. Having said that, any successful scenario must 
also explain the optical properties of the underlying stellar 
component. For this purpose, we examine the position of 
our simulated target galaxy in the CMR of observed ellipti- 
cal galaxies. The CMR is the simple and well-studied optical 
scaling relation for elliptical galaxies, and gives strong con- 
straints on the formation history of el liptical galaxies (e.g. 
lArimoto fc Yoshii lll987l:lGibson1ll997l) . On the other hand, 
due to our self-consistent treatments of chemodynamics, the 
optical properties of the simulation end-products can be de- 
rived with the population synthesis as described in Section 
12.3.21 and reflect the properties of stellar components, such 
as the age and the metallicity. Therefore, the CMR provides 
one of the best tests for our theoretical model of ellipti- 
cal galaxies. Fig. |5] shows the comparison of the simulated 
galaxies with the observed Coma cluster galaxies as well as 
NGC 4472 in theV - K sxiAU -V CMR. The data for the 
galaxie s in the Coma cluster are from lBower. Lucev. fc Ellis I 
119923). Since there is no difference between SOs and el- 
lipticals in the scaling relations which we discuss, we do 



' 1 ' J 
MoHpI R . 


• . /V . 


J 








1 1 h 

Model C / 


1 1 1 I 1 1 1 1 1 1 1 1 1 


MM 








,.•1 


: 





-28 -26 _-24 4 5 6 7 
Log Pg„3 (g cm ^) Log T (K) 



Figure 7. [Fe/H] vs density (left) and [Fc/H] vs temperature 
(right) distributions of gas particles within the projected radius 
of R =35 kpc for Model B (upper) and Model C (lower). The hor- 
izontal gray solid lines show [Fe/HJjf obtained by the spectrum 
fitting within the radius of ij = 35 kpc. The horizontal gray dot- 
ted lines correspond to mean [Fe/H] of stellar component within 
the radius of i? = 35 kpc. In the right panel, the vertical gray solid 
lines show the mean temperatures weighted by emission measure 
of the two temperature components (described by dotted lines) 
obtained by the spectrum fitting within the radius of i? = 35 kpc. 



not dis tinguish SOs from ellipticals. iBower. Lucev. fc Ellis I 
(|l992i) supplies the U — V and V — K colours which re- 
fer to an aperture size of 11 arcsec and the V band total 
magnitude derived from a combination of their data and 
the literature. We adopt the distance modulus of the Coma 
cluster of m — M = 34.7 m ag: the Virgo cluster distance 
modulus is m — M = 31.01 JCraham et al. Ill999ll and the 
relative distance modulus of Coma with respec t to Virgo is 
m - M = 3.69 iBower. Lucev. fc Ellis Ill992bl) . This gives 
a luminosity distance of 87.1 Mpc for Coma. We assume 
that the angular diameter distance equals the luminosity dis- 
tance, because the redshift of the Coma cluster [z ~ 0.023) 
is nearly zero cosmologically. Then the aperture size of 11 
arcsec at the distance of the Coma cluster corresponds to ~ 
5 kpc. Thus, we also measure the colours within an aperture 
diameter of 5 kpc for the simulation end-products. The to- 
tal luminosity of the simulation end-products is defined as 
the luminosity within the aperture radius of 50 kpc, which 
covers almost the whole simulated galaxy (Fig.|5J. 

We can see immediately that the colours of the result- 
ing stellar components of both Models B and C are too blue 
and inconsistent with the observational data. This incon- 
sistency can be traced to an excessive population of young 
and intermediate age stars. The strong feedback in Model 
C suppresses the formation of young stars, and mitigates 
this problem especially in the U — V CMR, which is sensi- 
tive to a population of young stars. However, Model C still 
has too blue a colour to reproduce the observed colours. 
Fig. 1^ shows that the histories of the star formation and 
SNe II and la rates for all the stars within r = 35 kpc. 
We can see that significant amounts of stars are continu- 
ously forming until z = 0. The SNe II and la rates are 
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Figure 8. Comparison of the simulated CMRs (filled circle and 
square for Models B and C, respectively) and that of the Coma 
cluster ellipticals (crosses). The dashed line shows the CMR fitted 
to the Coma Cluster galaxies. The large (small) open circle and 
square demonstrate the colour and magnitude for Models B and 
C, when the contribution from stars whose age is younger than 2 
(8) Gyr is ignored. Open star shows the position of NGC 4472. 



also too high. ICaperrallo. Evans, fc Turatto I il999l) pro- 
vide the observed SNe II and la rates of nearby ellipti- 
cal galaxies. Their SNe II rate gives the upper limit of 
0.02 X (Ho/75)^ SNu, and the rate of SNe la which they 
estimated is (0.018 ± 0.06) x (Ho/75)^ SNu. Here, Ho is 
the Hubble constant ( Ho = 70kms~^Mpc~^ is assumed 
in this paper) and SNu is the number of SNe per 100 yr 
and lO'°Ls 0, i.e. SNu = SN{imjr)-^{lQ^° Lg Q)-^ . In 
our models, the SNe II rate is significantly higher than 
SNe la rate, which obviously contradicts this observation. 
In Fig. IHl symbols with error-ba rs sho w the SNe la rates 
of lCaperrallo. Evans, fc Ttirattol ()l99flh . which is transfered 
from units of SNu to units used in Fig. |2l based on the B 
band luminosity for each model. Since the luminosities of 
both models are too bright, compared with the observed el- 
liptical galaxies (Fig. |HJ , these rates overestimate the obser- 
vational data. Nevertheless, both Models B and C obviously 
have too many SNe la. 

We also tried to fit the surface brig htness profile by the 
Sersic law (see eq. 11 of lKawata l2001a) . However, we could 
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Figure 9. Time variation of the star formation rate (upper) and 
the event rate of SNe II (thin lines) and SNe la (thick line) for 
Models B (left) and C (right). Solid circle (square) with error-bar 
in lower l eft (right) panel is taken from the ob servational SNe 
la rate bv ICaperrallo. Evans. &: Turatt'ol ^^9^. To show them 
clearly, those SNe la rates are plotted at t = 13 Gyr. 



not obtain any acceptable fit, because the stellar component 
is much too centrally concentrated. Such a central concentra- 
tion of stars is also induced by radiative cooling which leads 
to continuous star formation in the high density central re- 
gion. Instead, we measured the half-light radius as a radius 
where half the total luminosity is contained, when the total 
luminosity is defined as the luminosity within the aperture 
diameter of 100 kpc. Then, we get V band half-light radii 
of 4.2 and 3.4 kpc for Models B and C, respectively. These 
are too small to reproduce the observed Kormendy relation, 
the relation between the half-light radius and the luminosity. 
Therefore, our models reproduce neither the surface bright- 
ness profile nor the size of observed elliptical galaxies, owing 
to excessive central radiative cooling. 

Therefore, we conclude that although radiative cooling 
helps to explain the observed X-ray luminosity, temperature, 
and metallicity of elliptical galaxies, the resulting cooled gas 
also leads to the unavoidable overproduction of young and 
intermediate age stellar populations, at odds with the obser- 
vational constraints. Hence, the test for both X-ray and op- 
tical properties gives stronger constraints on the theoretical 
model. We claim that radiative cooling alone cannot explain 
both X-ray and optical properties, but another physical pro- 
cess to suppress the recent star formation induced by cooling 
is required. 

We found that even our extremely strong feedback as- 
sumed in Model C cannot stop radiative cooling which leads 
to the successive star formation. We compared the X-ray lu- 
minosity with the energy emitted by SNe II and SNe la in 
Model C, within r = 35 kpc at z=0. The X-ray luminosity in 



the 0.5-10 keV pass band is 3.3 x 10"^ erg 



The energy 



emitted from SNe is estimated using the mean SNe rate for 
the last 1 Gyr, and provides 4.4 x 10*'^ erg s~^. The energy of 
SNe seems larger than the X-ray luminosity. However, this 
X-ray luminosity is the luminosity in the 0.5-10 keV pass 
band. We measured the X-ray luminosity in the O.leV and 
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1 MeV pass band, and obtained 2.0 x 10''^ erg s^^. This is 
consistent with the total coohng rate, 1.1 x 10*^ erg s~^, 
of the gas particles within r — 35 kpc, which is calculated 
with the cooling rate used in our simulation code. Hence, 
the total cooling rate is much higher than the SNe feedback 
energy. Model C has a higher SNe rate than the observed el- 
liptical galaxies (Fig. |5J , and the X-ray luminosity of Model 
C is consistent with the observation (Fig.^J- Thus, this re- 
sult means that in the real elliptical galaxies SNe la do not 
generate enough energy to stop cooling. 

As expected, if the contribution of the young stars was 
ignored, the resulting colours match the observed CMR. 
Open circles and squares in Fig. |S| demonstrate that if 
the star formation in the last 2 Gyr did not occur, the 
CMR of both models would become close to observational 
data. If the stars formed in the last 8 Gyr did not ex- 
ist, the observed CMR would be roughly recovered by 
the increase in the colour and the decrease in the lumi- 
nosity. Therefore, one exotic solution is to "hide" these 
younger stars within a bottom-heavy initial mass func- 
tion (IMF) such that they cannot be scon today even 
if they did exist (e.g. Fabian, Nulscn, & CaizaresJ Il982t 
Ikathews fc Brighenti lll999D . 

Another (more plausible) possibility is that extra heat- 
ing sources, such as intermittent active galactic nuclei 
(AGN) activity, suppress star formation at low redshift. Be- 
fore suggesting this is the true solution though, we must 
re-examine the predicted X-ray properties of the simula- 
tion end-products after introducing these additional heat- 
ing sources; we will be pursuing this comparison in a future 
paper. 
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Figure 10. Radial dependences of Tx and [Fe/HJjf for Mod- 
els B (left panels) and C (right panels). Different symbols cor- 
respond to the values derived by the different methods, such as 
spectrum fitting (circles). X-ray luminosity (0.5—4 keV) weighted 
mean values (diamonds), p^T^^"^ weighted mean values for the gas 
with T > 12000 K (triangles), and p^T^^"^ weighted mean values 
for the gas with T > 10^ K (stars). In the upper panels, the 
error-bars show the two temperatures used in fitting the X-ray 
spectrum, and circles indicate the mean temperatures weighted 
by the emission measure of the two components. The horizontal 
error-bars show the radius range, and the vertical error-bars indi- 
cate the 90 % confidence range obtained by the spectrum fitting. 



3.3 Radial Dependences of X-ray Properties 

In Section 13.11 we focused on the mean properties in the 
central region. In this section, we discuss X-ray properties 
at different radii. First, we briefly mention the dependence 
of the results on the analysis method. Fig. llUl shows the pro- 
jected radial dependences of Tx and [Fe/H]x obtained by 
the following four different methods. Method 1 is spectrum 
fitting (plotted by circles and error-bars ), which is used in 
Section 13.11 method 2 takes the mean values weighted by 
the X-ray luminosity predicted by vmekal in the 0.5 — 4 keV 
energy range (diamonds); method 3 takes the mean values 
weighted by p^T^/^ for the gas with T > 12000 K (trian- 
gles); method 4 is the same as method 3 but only for the gas 
with T > 10^ K (stars). We can see that the spectrum fit- 
ting gives remarkably similar temperatures and abundances 
to method 4. Therefore, the results derived by the X-ray 
spectrum fitting reflect the properties of the high tempera- 
tures (r >~ 10^ K) gas. This is due to the X-ray pass band 
which is applied in fitting the spectrum (in this paper, 0.5—4 
keV, which is a typical ra nge used in the observational stud- 
ies). As shown in Fig. 1 of lMathiesen fc Evrardi J200ll) . the 
shape of the spectrum depends on the temperature of the X- 
ray emitting hot gas. As a result, the gas with T <~ 10® K 
hardly contributes to the X-ray spectrum at energies higher 
than 0.5 keV. For the same reason, method 3 provides dif- 
ferent values from the ones obtained by method 1. In partic- 
ular, the temperature of the central region is well bellow the 
ones obtained by method 1, due to the overestimate of the 
contribution from the low temperature (T <~ 10^ K) gas. 



Thus, if the simple p^T^^'^ weighted mean value is used to 
compare the simulation results with the X-ray observational 
data, the temperature limit for the X-ray emitting hot gas is 
crucial, and should be set a temperature as high as T = 10® 
K. 

Except in method 3, the temperature gradient has neg- 
ative slope. This is inconsistent with the observed temper- 
ature gradients in elliptical galaxies, which has a positive 
gradi ent especially at small radii fe.g. iFinogxienov fc .Tones I 
1200(31 . The observed positive gradient in temperature is con- 
sidered to be evidence for a cooling fiow. On the other hand, 
since the tcooi = tu line has positive slope in the density vs. 
temperature diagram (the right panels of Fig. (SJ, including 
radiative cooling leads to a negative temperature gradient, 
if the gas density profile also has a negative gradient, which 
is naturally expected. Thus, to explain the observed posi- 
tive temperature gradient, i.e. the low temperature (T <~ 1 
keV) in the central region, significant amounts of the gas has 
to stay in the region of tcooi < Ih in the density vs. temper- 
ature diagram for a long time, ~ Ih- One possible physical 
process to realise this condition is again heating sources to 
balance radiative cooling, although it requires a fine-tuned 
physical process. In fact, the stronger feedback in Model C 
leads to the shallower temperature gradient than Model B 
(see also Fig. (SJ, although Model C still fails to reproduce 
the positive gradient. 

In [Fe/H]x, irrespective of the method, there is a 
clear negative gradient, which is consistent with the ob- 
servational trends of som e bright elliptical galaxies (e.g. 
iFinoguenov fc Jones l2000h . This is because stars, which are 



14 D. Kawata and B.K. Gibson 



the source of the iron, stay within the central region cor- 
responding to the inner most bin. This explains why the 
smaller aperture gives a higher [Fe/H]x in Fig. El 

In [Fe/H]x, method 1 gives similar results to method 
4. Particularly in the central bin, method 3 gives a higher 
metallicity than method 1, due to the contribution from the 
cooler {T < 10® K) gas which is in the central region (Fig.|^ 
and enriched by stars preferentially. Surprisingly, method 2 
gives higher metallicity than method 1. We found that this 
is because the amount of line emission increases with metal- 
licity, and method 2 is biased towards higher metallicity gas. 
In fact, if the continuum X-ray luminosity is used, and line 
emissions are ignored, the luminosity weighted mean val- 
ues agree with the values obtained by the spectrum fitting. 
Therefore, we should be careful when we compare the results 
of theoretical models with the X-ray observations. 

Top panels in Fig. 1111 show the 0.5-4 keV X-ray lu- 
minosity profile. We fit these profiles with the foUwoing j3 
model, which is used to fit the observed X-ray luminosity 
profiles (e.g. iFor man. Jo nes, fc Tucker I llQSSt iMatsushita I 
ll997l:lMatsushita et al. Ill998h . 

Ex {R) = Sx,o [l + {R/Rco..f] . (3) 

In fitting, we used the data at radii less than 100 kpc, be- 
cause the X-ray luminosity drops rapidly at radii further 
than 100 kpc. Thin lines in Fig. 1111 show the best fit pro- 
file, and the parameters are summarised in Table |1] The 
X-ray luminosity profiles for all the Models are well de- 
scribed with the /3 model. The typical values in the ob- 
served elliptical galaxies are /3 ~ 0.5 and -Rcore ^ 1 kpc 
l|Forman. .lones. fc Tuckerl I198RI : iMatsushital Il99/t l. The 
values of /3 obtained for our models are similar to the ob- 
served values. It is worth noting that the profile for Model A 
is significantly steeper than those of Models B and C. Thus, 
radiative cooling leads to a shallower X-ray luminosity pro- 
files, i.e., lower f3. This is also a consequence of the decrease 
in the central hot gas density with radiative cooling (Fig. 

On the other hand, the core radius is too large to repro- 
duce those for the observed elliptical galaxies, irrespective 
of models. This is also a problem to solve for the current nu- 
merical simulation model. In model A, the core comes from 
the central region where the density profile is fiat due to 
the thermal pressure of the adiabatic gas (Fig. |5J . On the 
other hand, in Models B and C, the large core radii might 
be caused by a poor resolution of the simulation, because in 
the central region the number of the hot gas is small due to 
radiative cooling. 

The X-ray luminosity profile is used to estimate the 
total mass, Mtot(»'), within a three-dimensional radius, r, 
combined with the temperature and metallicity profiles (e.g . 
iForman. Jones, fc Tucker I Il985l : IMatsushita et al. 1 119981 . 
Assuming that the X-ray emitting hot gas is in a hydro- 
static equilibrium and has a spherically symmetric distribu- 
tion, Mtot,hyd(»') is written by 



Mtot,hyd(r) = -- 



fcTg(r)r /rflnpg(r) dlnTg(r) 



dlnr 



+ 



dlnr 



(4) 



where G, nip, and k are respectively the constant of gravita- 
tion, the proton mass, and the Boltzman constant; Tg and Pg 
are the temperature and density of the X-ray emitting hot 
gas; jj. = 0.6 is the mean molecular weight of the gas. Since it 
is easy to analyse Mtot (r) from the output of numerical sim- 



ulations, it is interesting to compare the mass profile which 
is derived using equation @, Mtot,hyd(r), with that derived 
from the simulation outputs directly, i.e. the summation of 
the mass of particles within the radius, Mtot,sum(»'). To use 
equation Q, we need the hot gas density and temperature 
profile. We derived the hot gas density using the relation of 

[ex(r)/;n,ekai(r(r),Z(r))]i/^ 



(5) 



Here, ex is the X-ray emission per volume, which 
is calculated from the f ollowing Abel integration (e.g. 
iBinnev" fc Merrifield Il998^ of the surface luminosity profile, 
Ex(-R), of equation 



ex{r) 



dEx(-R) dR 



dR /R2 



(6) 



We numerically integrate this equation and obtain exij')- 
The expected X-ray emission, /mokai (T,Z), is a function 
of temperature and metallicity, and calculated from the 
XSPEC vmekal model in the same X-ray pass band as 
Ejf (0.5-4 keV). Here, we use the metallicity, Z, which is 
not [Fe/H], but the total abundance of heavy elements, 
and assume the meteoric solar abundance pattern. Next, 
the three dimensional temperature and metallicity profiles 
are required. In observational studies, the three dimen- 
sional temperature and metallicity profiles are often as- 
sumed to equal the projected two dimensional profiles. Re- 
cently, the deprojected three dimensional profiles are also 
derived, assuming a s pherical symmetry, and used to esti- 
mate Mtot,hyd(?') (e.g. IMatsushita et aLll2002l) . Therefore, 
we consider the two cases of the projected profile and that 
derived from the three dimensional distribution of particles 
(deprojected profile). The thick lines in the panels of the 
second and third rows of Fig. Illl show the temperature and 
metallicity profiles of the X-ray emitting hot gas for all the 
models. Here, we use method 4 to derive the profiles, be- 
cause methods 1 and 4 give similar X-ray observed values, 
as shown above. Both temperature and metallicity profiles 
are fitted with a power-law function, which is written by 
Tg(J?) = rg,o[l -f (R/Rt )]""^ in the case of the projecte d 
temperature profile (e.g. IForman. Jones, fc Tuckerl Il985^ . 
The best fit profiles are presented as thin lines in Fig. 1111 
and the parameters are summarised in Table (31 In fitting, 
we used the data at radii less than 100 kpc, which is the 
same range as that used in fitting of the luminosity profiles. 
Only for the deprojected temperature profiles of Model C, 
we exclude the inner most bin whose value is obviously in- 
consistent with the assumed profile. In Model A, the metal- 
licity is set to be zero. The fitting of the metallicity profiles 
in Models B and C is poor, because the profiles are not 
smooth. However, the metallicity profiles are less sensitive 
to Mtot,hyd(?')- Finally, the bottom panels of Fig. 1111 shows 
the derived Mtot,hyd(?') as well as Mtot, sum (r) which is the 
summation of the mass of all the components, i.e. DM, gas, 
and stars, within the three dimensional radius, r, in the sim- 
ulation end-products. In Model A, Mtot,hyd(7') is consistent 
with Mtot,sum(''), which is because the X-ray emitting hot 
gas in Model A is in a hydrostatic equilibrium. In Models 
B and C, around the radius of 50 kpc, Mtot,hyd('~) matches 
with Mtot,aum(T'). It meaus that even if radiative cooling is 
included, the status of the X-ray emitting hot gas is close 
to a hydrostatic equilibrium in this region. It is not surpris- 
ing that Mtot,hyd(?') does not reproduce Mtot, sun (r) at larger 



Multi-wavelength Cosmological Simulations of Elliptical Galaxies 15 



radii (r > 100 kpc), because we exclude those data in fitting. 
On the other hand, at small radii Mtot,hyd(r) is systemat- 
ically smaller than Mtot,sum(r), because the hot gas is not 
in a hydrostatic equilibrium in this region, due to radiative 
coohng which continuously converts the hot gas to the cold 
gas. Since our simulations do not reproduce all the observed 
properties, we cannot judge whether or not this happens in 
the mass estimation from the observational data. It is also 
worth noting that Mtot,hyd('') from the deprojected profiles 
(dashed lines in Fig. IIH gives higher mass and reproduces 
Mtot,sum('') better than Mtot,hyd(r) from the projected pro- 
files (solid lines in Fig. IIH in Models B and C. This is be- 
cause the projected temperature underestimates the depro- 
jected temperature, when there is a significant temperature 
gradient (Fig. Since the temperature is constant inde- 
pendent of radius, in Model A, the projected temperature 
profile is consistent with the deprojected temperature pro- 
file, 

3.4 Time Variation of X-ray Properties 

In this paper, we focus on the properties of the simulated 
galaxies at 2: = 0. The advantage of numerical simulation 
is to be able to study the evolution of galaxies, while ob- 
servational data provide a snapshot of galaxies. Here, we 
briefly discuss the evolution history of the simulation end- 
products. Fig. ll2l shows the X-ray properties and the hot gas 
and baryon mass of the simulated galaxies within ii = 35 
kpc as a function of time. The X-ray luminosities are mea- 
sured in the 0.5-10 keV pass band, and temperatures and 
iron abundances are measured by method 4 in Section [3.31 
Since the number of output flies in the simulations are lim- 
ited, the lines are not smooth. In the X-ray luminosity, all 
the models have the relatively constant evolution, except 
for temporary increases. These temporary features are asso- 
ciated with the accretion of small bound objects, i.e., minor 
mergers. In fact, we confirmed that a small bound object 
passed through the central region of the target galaxy, and 
the gas of the objects is striped, when the X-ray luminosity 
increases temporarily in Model A. It is also seen in Fig. 1121 
that at that time the baryon and hot gas mass increases, 
and the temperature also increases. Although the X-ray lu- 
minosities in Models B and C are not smooth, the X-ray 
luminosity is always higher in the order of Models A, C, 
and B, as seen in Section 13.1.11 On the other hand, there 
is significant fluctuation in the temperatures for Models B 
and C, and a trend seen in Section r3.1.1l i.e. Model B has a 
higher temperature than Model A, represents a condition in 
the last 1 Gyr, when the systems are relatively stable. The 
mean trend of the iron abundance is constant. However, the 
iron abundance also changes, and gets higher values at the 
same time as the temperature. We consider that this is due 
to the SNe feedback followed by the minor mergers, which 
heats and enriches the gas. The baryon mass within this 
aperture increases continuously in Models B and C, because 
the gas turns into stars and the fresh gas accretes from the 
outer region. Since the hot high pressure gas stays in Model 
A, the baryon mass is relatively constant. 

The numerical resolution is always a serious issue for 
the galaxy formation simulation, although we have achieved 
a high resolution at a similar level to the recent other cos- 
mological simulations which include cooling and star forma- 



tion and follow the evolution until z = (e.g. iMeza et al. I 
2003). Since we focus on the mean properties in the scale 
larger than our spatial resolution, we believe that our con- 
clusion is not heavily affected by the resolution. To examine 
this, we carried out a ~ 3.375 times (by mass) higher res- 
olution simulation with the same parameters as Model B. 
In an initial condition for the higher resolution model, we 
have assigned high-resolution particles only in the region 
surrounding the target galaxies (roughly within the radius 
of 4x Tvir at z= 0), to reduce the computational time. There- 
fore, the random seed for the small density perturbation is 
different from Model B, which leads to a different history of 
minor mergers. The mass and softening length of each gas 
(DM) particle in the high-resolution region are 1.74 x 10^ 
(1.17 X 10**) Mq and 1.51 (2.86) kpc, respectively. We stop 
the simulation at z ~ 0.25, when low-resolution particles 
enter the virial radius of the target galaxy. The gray lines in 
Fig. 1121 show the evolution of the X-ray properties and the 
hot gas and baryon mass of the simulated galaxies within 
-R = 35 kpc for the high-resolution model. Although detailed 
features are different, the mean values are roughly consistent 
with those of Model B. Hence, we conclude that the results 
presented in this paper are not sensitive to the numerical 
resolution. However, even in this high-resolution model the 
resolution is still much larger than the scales of the various 
physical processes, such as the star formation and the SNe 
feedback. The effect of such small scale physics on the large 
scale properties is little known, and it depends on how to 
model those processes. Therefore, the results shown in this 
paper might be affected by the modeling we adopt. Needless 
to say, this modeling of the small scale physics is the most 
crucial issue for studies of galaxy formation, which requires 
further investigations. 



3.5 Abundance Ratios of the X-ray Emitting Hot 
Gas 

Since we follow the evolution of the abundances of various 
heavy elements, we can analyse the abundance ratios of dif- 
ferent elements. Especially, the ratios of a elements and Fe 
provide useful information about the chemical enrich history, 
because a elements, such as O and Si, are mainly produced 
by SNe II, while Fe is mainly produced by SNe la. Table |S] 
shows the abundance ratio of O and Si with respect to Fe 
(0/Fe and Si/Fe) for the X-ray emitting hot gas, the cold 
gas, and the stellar component within the radius of R =35 
kpc. For the hot gas, the values obtained by method 1 and 
method 4 are shown. For the cold gas and stars, the mean 
values weighted by the mass are presented. We found that 
Si/H and O/H measured by method 1 are a little different 
from the values measured by method 4, which causes dif- 
ferences in 0/Fe and Si/Fe between methods 1 and 4. We 
consider that this is because of a difficulty of the X-ray spec- 
trum fitting due to line blending. Nevertheless, both meth- 
ods give less than solar O/Fe and Si/Fe. This is consistent 
with the recent XMM-Newton RGS observations of O/Fe 
and Si/Fe in the central region of elliptical galaxies (e.g. 
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Figure 11. Top row: the 0.5-4 keV X-ray radial luminosity profiles (thiclf lines) and the best-fit beta profile (thin lines). 2nd row: 
projected (solid lines) and deprojected (dashed lines) temperature profiles (thick lines) and the best fit profiles (thin lines). 3rd row: the 
projected (solid lines) and deprojected (dashed lines) metallicity profiles (thick lines) and the best fit profiles (thin lines). Bottom row: 
Mtot,sum(r') (thick lines) and Mtot,hyd('') (thin lines). The solid and dashed lines corresponds to the profiles calculated firom projected 
and deprojected temperature and metallicity profiles, respectively. In model A, the solid and dashed lines are almost overlapped, because 
they provides similar values. 



Table 4. The results of the fitting of the luminosity, temperature, and metallicity profiles. The values in round brackets denote the 
results of the fitting of the deprojected profiles. 





log5o 


H 


-f^core 


logTx,o 




-RT(fT) 


logZo 


az 


Rz{rz) 


Name 


(erg s~^ pc~^) 




(kpc) 


(keV) 




(kpc) 


{%) 




(kpc) 


A 


33.6 


0.86 


23.8 


-0.331 (-0.328) 


0.15 (0.15) 


128 (141) 








B 


31.5 


0.35 


16.8 


0.250 (0.327) 


0.63 (0.88) 


26.8 (47.2) 


-1.14 (-1.07) 


0.67 (1.80) 


26.2 (98.7) 


C 


32.5 


0.51 


20.4 


0.097 (0.270) 


3.88 (0.88) 


399 (39.9) 


-0.38 (0.19) 


0.83 (0.98) 


2.36 (1.32) 
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Figure 12. Time variation of the X-ray luminosity (upper left), temperature (lower left), iron abundance (upper right), and the hot 
(thin lines) and baryon (thick line) mass (lower right) within R = 35 kpc. The values obtained in Model A/B/C are presented as 
dotted/solid/dashed lines. In the lower right panel, since in Model A almost all the baryon component within this aperture is the hot 
gas, the thick line overlap with the thin line. The gray line shows the results of a higher resolution simulation with the same cooling and 
feedback model as Model B. 



O/Fe=0.4 solar in N GC 4636^ O/Fe=0.37 and Si/Fe=1.0 4 
solar in NGC 5044: IXu etlTI 120021: iTamura et al. 1120031) . 

Table |S| shows the mass fraction of Fe, O, and Si ejected 
from SNe II, SNe la, and the IM stars within the last 1 Gyr 
for Models B and C at i? < 35 kpc. Table El also presents 
O/Fe and Si/Fe of the total ejected gas. It appears that 
SNe la are a dominant source of Fe, which induces the low 
O/Fe and Si/Fe. Higher Si/Fe than O/Fe is also explained 
by a significant contribution of SNe la which produce a fair 
amount of Si, compared to O. However, note that the rates 
of SNe II and SNe la are too high in both Models B and C, as 
discussed in Section r3.2l (Fie.l^. Especially, the SNe II rates 
in the observed elliptical galaxies are very small, and thus 
our high SNe II rate overestimates O and Si production. In 
addition. Table |K| shows that O/Fe and Si/Fe for the stellar 
component are solar or less, while the observed bright el- 
liptical galaxies have over solar abundance ratios (Q/Fe~2). 

^ We assumed that lXu et al. I j20o2) used the photospheric solar 
abundances, and converted their results to the values with respect 
to the meteoric solar abundances. 



This is also because stars continuously form from the cold 
gas which has low O/Fe and Si/Fe (Table 0. This low a/Fe 
for the stellar component underestimates O/Fe and Si/Fe for 
the mass-loss from the IM stars, although the contribution 
of the mass-loss of IM stars to the metal enrichment is small, 
compared with SNe. Hence, strictly speaking, the chemical 
enrichment history of our simulation is inconsistent with the 
observation. Nonetheless, broadly speaking, our simulation 
demonstrates that the observed low Si/Fe and O/Fe in the 
hot ISM of elliptical galaxies are explained by the dominant 
contribution of SNe la to the chemical enrichment of the 
X-ray emitting hot gas. 



4 DISCUSSION AND CONCLUSIONS 

Our cosmological chemodynamical simulations make it pos- 
sible to undertake quantitative comparisons between the 
theoretical models and the observational data in both the 
X-ray and optical regime with minimal assumptions. This 
is the first attempt to explain both the X-ray and optical 
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Table 5. Oxygen and Silicate abundances with respect to iron for X-ray emitting hot 
stellar component at R <35 kpc. 



'as, cold gas, and 



Name 
B 

C 



hot gas 



method 1 
O/Fe Si/Fe 
(solar) (solar) 
0.21 0.22 
0.91 0.64 



method 4 
O/Fe Si/Fe 
(solar) (solar) 
0.20 0.51 
0.60 0.82 



cold gas 
(T< 10^ K) 
O/Fe Si/Fe 
(solar) (solar) 
0.63 0.91 
0.51 0.80 



stars 



O/Fe 
(solar) 
1.02 
0.89 



Si/Fe 
(solar) 
1.13 
1.07 



Table 6. Ejected heavy element mass fraction from SNe II, SNe la, and IM stars and the abundance ratios 
of the total ejected gas. The values denote a mean values within the last 1 Gyr from 2 = at i? < 35 kpc. 



SNe II 

fcj.Fc" foj.O 

B 0.29 0.74 
C 0.19 0.62 



fcj.Si 

0.52 
0.40 



fcj.Fc 

0.50 
0.63 



SNe la 

fcj.O 
0.02 
0.04 



fcj.Si 
0.25 
0.36 



fcj.Fo 

0.21 
0.19 



IM stars 

0.24 
0.34 



fcj.Si 

0.23 
0.24 



total 
O/Fe Si/Fe 
0.62 0.90 
0.48 0.77 



fn 



cj.Fo ■ 



M, 



cj.Fc.SNc 



Il/Mej, 



properties of observed elliptical galaxies by self-consistent 
cosmological simulations. First, we found that radiative 
cooling is required to explain the observed X-ray lumi- 
nosity, temperature, and metallicity of elliptical galaxies. 
Comparison between models ignoring and including radia- 
tive cooling clarifies that radiative cooling ensures that the 
hot dense gas turns into the cold (i.e. non X-ray emit- 
ting) gas, and keeps X-ray emitting gas at high temper- 
ature and low density (Fig. (SJ. As a result, including ra- 
diative cooling leads to lower X-ray luminosity and higher 
X-ray weighted temperature, and provides a better agree- 
ment with the observed Lx — Tx relation (Fig.|lJ. Although 
this effect of c ooling has already boon sho wn by previ- 
ous studies (e.g. lPearce et al. 2 000: Muanwong"et al. 120021: 
iDave. Katz. fc Weinberg 1 12002|) . our study newly confirms 
that radiative cooling works in the same way on galactic 
scales using simulations with higher physical resolution. In 
addition, we find that stronger SNe feedback leads to higher 
X-ray luminosity and lower temperature, because heating 
by the SNe feedback makes the cooling time of the hot 
gas longer, and thus the gas with higher density and lower 
temperature is allowed to stay as X-ray emitting hot gas 
(Fig. |SJ . Radiative cooling also ensures that the hot gaseous 
halo has not been enriched efhciently. Stars preferentially 
enrich the gas in the central region, where cooling is effi- 
cient. The enriched gas can then cool easily and be incor- 
porated into future generation of stars. In fact, we found 
that a large fraction of iron ejected from stars is locked into 
stars. This effect of cooling has been alrea dy suggested by 
iFuiita. Fukumoto. fc Okoshi I (Il996l Il997l) . ahhough they 
did not consider any star formation induced by the cooled 
gas. Our more self-consistent numerical simulation confirms 
their scenario. 

To summarise the X-ray study, our radiative cooling 
models succeed in reproducing the X-ray luminosity, temper- 
ature, and metallicity of observed elliptical galaxies. How- 
ever, we also found that the resulting cooled gas leads to 
unavoidable overproduction of young and intermediate age 
stellar populations, at odds with the optical observational 
constraints. We examined the position of our simulated 
galaxy in the observed Coma cluster CMR, and found that 
the colours of the resulting stellar components are incon- 



sistent with the observational data (being too blue. Fig. (HJ. 
Therefore, this study demonstrates that the cross check over 
both X-ray and optical properties gives stronger constraints 
on the theoretical models, and is essential for any success- 
ful scenario of elliptical galaxy formation and evolution. For 
example, our results throw doubt on the recent claim that 
only radiative cooling is required to explain steeper slope of 
the Lx — Tx relations of groups and elliptical galaxies. 

Consequently, our model fails to reproduce both X-ray 
and optical properties of the observed elliptical galaxies. 
Nevertheless, our study provides some hints to explain those 
observed properties. The biggest problem of the present 
model is too many recent (z<l) star formation (Fig. (HJ. 
Therefore, one exotic (if somewhat ad hoc) solution is to 
"hide" these younger stars within a bottom-heavy initial 
mass function (IMF) suc h that they cannot be seen today 
even if they did exist (e.g. [ Fabian. Nulsen. fc Caiz arcs 1983; 
iMathews fc Brighenti Ill999t) . A more plausible possibility 
is additional physical process which suppresses the cooling 
and star formation after the system has formed. Although 
we assume the extremely strong SNe feedback in Model C, 
it is not enough to suppress the star formation Hence, the 
SNe la feedback is unlikely to be a candidate of this phys- 
ical process, although it should contribute in some degree. 
We expect that extra heating by intermittent AGN activity, 
which is not included in our numerical model, is a major 
mechanism to suppr ess the star formation at low red shift 
jBrighenti fc Mathews l2002l: iBriiggen fc Kaiser l2002h . Re- 
cent observations suggest that all elliptical galaxies have a 
central black hole. Therefore, once the cold gas falls into 
the central part of galaxy, the AGN is expected to be ac- 
tive and heat up the surrounding gas, and might be able to 
blow out the gas from the system via a radio jet. Once the 
gas accretion is stopped by AGN heating, the AGN itself 



* Note that this might come from our poor modeling 
of the SNe feedback. The energy feedback from SNe is 
one of the most difficult processes to model in galaxy 
formation simulations. Various SNe feedback models hav e 
been suggeste d (e.g. |NaAa£ni.lL-WhiteJ E,99j; |Y cpes et al 
1997; Hultma n fc Pharasvn lll99Sl: iThacker fc Couch nian 200* 



El 



9 



Springel fc Her"nQuis^l20o3v Unfortunately they neglect SNe la. 
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loses fuel and goes back to a quiescent state, which allows 
the gas to cool again. This self-regulation cycle might be 
able to suppress the star formation significantly for a long 
time. In fact, recent high resolution X-ray images taken by 
Chandra reveal that the central hot gas in elliptical galaxies 
is not smoothly distributed, but is c avitated on scales com- 
parable to the radio emiss i on fe.g. iBohringer et aTlllQQd 
iFinoeuenov fc Jones lf200ll: I Jones et al. 1 1200^. The X-rav 
holes are also seen on the galaxy clusters scale (< 50 kpc), 
and are often coincident with extended radio lobes (e.g. 
i McNamara et al. I I2OO0I: lAUen et aTl I2OOII : iBlanton et aU 
l20(UrEeint^^dnl20n2^ . Those features are considered to 
be a relic of the significant infiuences of AGN activity on 
the hot gas of elliptica l galaxies and intra cluster medium 
JChurazov et al. ll2nnl. 

This AGN heating mechanism might also help to 
solve the problem in the observed [ Mg / Fel -magnitude 
relation of elliptical galaxies (e.g . jM^rteuccJ_^94 
[Mgj jteucci. Ponzone. fc Gibson I Il99?j : iKawata fc Gibson I 
l2003l) . To explain the [Mg/Fe]-magnitude relation, larger 
galaxies are required to stop star formation at an earlier 
epoch. Since the recent observations claim that the black 
hole mass increa ses with the velocity dispe rsions of their 
host bulges (e.g. iMerritt fc Ferrarese I I2OO1L and reference 
therein), the heating effect by AGN is expected to increase 
with the stellar mass of elliptical galaxies, and star forma- 
tion in the larger galaxies are suppressed more strongly once 
the AGN acts. We are planing to construct self-consistent 
models of this mechanism, and examine how such a self- 
regulation mechanism induced by AGN affects the X-ray 
and optical properties for elliptical galaxies with different 
luminosities. 
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